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Abstract. We show that an enslaved phase-separation front moving with diffusive 
speeds U = C/ VT can leave alternating domains of increasing size in their wake. We 
find the size and spacing of these domains is identical to Liesegang patterns. For equal 
composition of the components we are able to predict the exact form of the pattern 
analytically. We also show that there is a critical value for C below which only two 
domains are formed. Our analytical predictions are verified by numerical simulations 
using a lattice Boltzmann method. 



1 Introduction to Liesegang Patterns 

The formation of highly ordered patterns in naturally occurring biological, chemical, and 
mineralogical systems has long been a subject of intense interest. The study of such pat- 
tern formation can sometimes allow deep insight into their underlying natural phenom- 
ena. In a previous paper we analyzed the dynamics of pattern formation behind a one- 
dimensional, slow moving (enslaved), phase-separation front. Our analysis concerned 
fronts moving with constant speed, and the pattern formed was a series of alternating 
bands of regular width and spacing [4]. We show here that a front moving with diffusive 
speed will form a more complex Liesegang pattern. 

It was with the motivation of understanding pattern formation in simple systems that, 
just over one hundred years ago, R. E. Liesegang observed a highly ordered pattern of 
concentric rings precipitating around a drop of silver nitrate on a glass slide with a thin 
gel coating containing potassium dichromate [8]. These concentric rings are now known 
as Liesegang rings. The radial width and spacing of the rings increases with increas- 
ing distance from the center. Rings close to the center are narrow and tightly packed. 
Rings far from the center are wide and far apart. The pattern forms from the center out- 
wards, and is stationary once visible. Several alternatives to silver nitrate and potassium 
dichromate in the production of Liesegang rings have been used in the literature. In gen- 
eral, some electrolytes A and B combine to form an insoluble precipitate D which then 
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produces Liesegang patterns. Most recent publications use linear Liesegang patterns of 
bands and gaps which are produced by adding the A electrolyte to a test-tube containing 
the B electrolyte suspended in gel [9]. 

To characterize Liesegang bands or rings, they are typically numbered from the first 
formed to the last formed, the nth appearing at time t n at position x n with a width w n . 
Repeated careful measurements of Liesegang patterns revealed that discrete, defect-free 
bands could be characterized by a set of empirical laws [6]. They are 

Time Law x n cx \ft~ n , (1.1) 
Spacing Law x n +i/x n = l+p, (1.2) 
Width Law x n <xw n . (1.3) 

The time law relates the position of the nth band with the time of its appearance. The 
location of subsequent bands is given by the spacing law, where p > is the spacing coeffi- 
cient. The width law states that the band width is proportional to the position of the band, 
which is a natural result of the spacing law with the assumptions of mass conservation 
and uniform concentration of precipitate bands [7]. These laws are only considered valid 
for large n. Much attention has been paid to the phenomenological Matalon-Packter law 
p=F(£>o) + G(£>o) / -^0/ which relates p to the initial concentrations Aq and Bo of the A and 
B electrolytes [2]. F(£>o) and G(£>o) are known to be decreasing functions of Bo, though 
not much else about them is known. Antal et al. [2] show that the Matalon-Packter law 
can be derived in limiting cases from more general expressions. 

Soon after their characterization, study into the nature and cause of Liesegang pat- 
terns was prolific. For instance, a paper by Stern in 1954 makes mention of more than 
six hundred papers having been published on the subject by that time [10]. There were 
many early attempts to develop a comprehensive model of Liesegang pattern formation. 
However, it proved difficult to account for the wide variety and complexity of possible 
patterns, and many early models were eliminated by additional experiments. The com- 
plexity of Liesegang pattern formation thwarted theoretical understanding, and progress 
slowed. More than a century later there is currently still no generally accepted compre- 
hensive mechanism for Liesegang pattern formation. 

Current dominating theories of Liesegang pattern formation can be categorized as 
either an ion-product supersaturation theory where electrolytes combine directly into the 
precipitate (A+B — > D), or a nucleation and growth theory where one or more intermediate 
compounds form before final precipitation (A + B — > C — » D). A brief discussion of these 
models is given in the recent paper by Jahnke and Kantelhardt [7]. 

For either of these theories, precipitation occurs behind a reaction front formed by 
the A electrolyte diffusing into the B electrolyte. Since the initial concentration (Aq) of 
A in the drop is typically an order of magnitude higher than the concentration (Bo) of 
B in the gel, the time dependent concentration profile A(x,t) of A electrolyte in the gel 
forms a reaction front which resembles the familiar heat-diffusion profile moving with 
speed u r cxt^ 1 ^ 2 [2]. Coincident with, or trailing behind this reaction front is a precipita- 
tion front, moving at a speed Up oc t~ 1//2 , which produces the banded pattern. The precise 
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relationship of the precipitation front to the reaction front depends on the specific chem- 
ical and physical mechanism of Liesegang pattern formation which are still not entirely 
understood. 

The recent trend for publications on Liesegang patterns shows a re-surging inter- 
est. A search of the ISI Web of Science Internet database for publications on the topic 
"Liesegang" shows an increase in papers since a fall in the early second half of the 20th 
century. We found only 8 papers from 1970 to 1979, 35 from 1980 to 1989, 107 from 1990 
to 1999, and 178 since 2000. Much of the recent research into Liesegang pattern forma- 
tion has focused on the moving front by employing a variety of numerical techniques 
to simulate models of precipitation fronts. Some examples are the reaction-diffusion cel- 
lular automata simulations by Chopard et al. [3], direct simulation of a model-B system 
with a chemical reaction like source term by Antal et al. [1] and Racz [9], the discrete 
stochastic simulation which used random walkers to model the diffusion front by Izsak 
and Lagzi [5], as well as the lattice gas simulation by Jahnke and Kantelhardt [7]. 

In this paper we show that patterns identical to Liesegang patterns can be formed in 
a much simpler physical system. We consider a binary mixture that can phase-separate if 
a control parameter crosses a critical value. In this system an enslaved phase-separation 
moves at a speed of u = cf~ 1//2 . An example of such a system would be a binary mixture 
that is cooled below its critical temperature from one end. This situation is somewhat 
similar to that of the electrolytes, where material is formed at the front and will sub- 
sequently phase separate. However the details are quite different, most notably in the 
nucleation conditions. 

It should be noted that this model is somewhat similar to the Model-B precipitation 
front proposed by Antal et al. [1], but there are several key differences. Their model 
has a moving, Gaussian shaped source of ^4-type material designed to mimic the prod- 
uct of a chemical reaction front. Their source moves through a region which is phase- 
separated into an equilibrium £>-type. When the concentration of „4-type material in a 
given area reaches the spinodal value, an „4-type domain nucleates and depletes the sur- 
rounding region of its ^4-type material. The source moves on, leaving stable domains. 
The speed, width, and concentration of their source are free parameters of their model. 
Our model has an abrupt control parameter front which induces phase-separation as it 
passes through a mixed material, the mechanism for nucleation of new domains is quite 
different. Notably switching does not occur at the spinodal value for reasons explained 
in [4]. Also the analysis of Antal et al. is numerical in nature. Our model has only one 
parameter and can be solved analytically. 

To allow for an analytical treatment we make the simplifying assumption of an abrupt 
front. We expect that the results will be qualitatively similar to those for an extended 
control parameter front. In practice, an abrupt front could be experimentally achieved by 
immersing a thin sample with a prescribed speed u <x f -1 ' 2 into a temperature bath. 

The key result is that we are able to analytically determine the resulting patterns and 
we show that they obey the Liesegang laws given by Equations {1.2) and \\.3) . We 
derive an analytical expression for the spacing coefficient p in terms of the free parame- 
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ters of this model. We verify these theoretical predictions by direct numerical simulations 
using a lattice Boltzmann method. 



2 A Model for Liesegang Patterns Formed by Enslaved Phase 
Separation Fronts 

We consider two materials, an ^4-type and a £>-type, in an incompressible mixture such 
that the total density p = p^(x,t) +pg(x,t) remains constant. The relevant variable is then 
the relative concentration of A to £>-type material defined as 



<p{x,t) 



p A {x,t)-p B {x,t) 



(2.1) 



From here on the time and position dependence of the concentration will be implied. For 
simplicity we assume that the two materials have the mixing free energy described by a 
(p 4 law 



dx 



a(x.t) 9 b(x,t) d K(x,t) ,_ 



(2.2) 



The time and position dependence of the control parameters (a, b and k) are such that they 
constitute a spatially abrupt transition from the mixing region to the phase-separating 
region of the phase diagram. For example: 



a(x,t)=a s + (a M -a s )© 



J udt+xo 



(2.3) 



where is the Heaviside step function, and the transition takes the form of a front mov- 
ing with velocity u(t). The free energy of Eq. \2.2\ has a single minimum for a > 0, result- 
ing in material mixing. When a < there are two minimums, resulting in the separating 
of material. The control parameters for the mixing and separating regions are denoted by 
subscripts M and S respectively. The other parameters are defined similarly. 

Since we assume incompressibility and a one-dimensional geometry of the system, 
hydrodynamics can be neglected here, and the dynamics is therefore purely diffusive 



d t (p = V[m(x,t)Vp{x,t)\ . 
The chemical potential is derived from the free energy as, 

p{x,t) = — =a(x,t)(p+b(x,t)(p 3 -K{x,t)V 2 (p, 



(2.4) 



(2.5) 



and the diffusive mobility m is one of the control parameters which can vary across the 
front similarly to (|2.3[) . Aside from time and space dependence of the control parameters, 
this is the familiar Model-B. 
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This model for a moving front has many parameters which can affect the dynamics of 
phase separation, however not all of these parameters are independent. As we elaborate 
on in our previous paper, introduction of appropriate time, space, and concentration 
scaling can non-dimensionalize the equations of motion [4]. For the remainder of this 
paper we will work entirely in the non-dimensional scales 

T=* X-f ,*=* (2.6) 

Here t sp = 4K/ma 2 and A sp = 2tt\/— 2k/ a are the characteristic time and length scales of 
spinodal decomposition, and (p eq =\ / —a/b is the positive equilibrium concentration of the 
phase-separated material. Non-dimensional quantities will be denoted by capital letters. 

Nondimensionalization reduces the free parameters to an independent set of four 
parameters: U = u/u sp is the speed of the front scaled by the natural speed of spinodal 
decomposition u sp =A sp / 't sp , M=niM/ms is the ratio of the diffusive mobility ahead of the 
front to behind the front, A — —cim I « s is the depth of the quench into the unstable region 
of the phase diagram, and <E>,>, = <pj n / (p eq is the non-dimensional initial concentration of 
the mixed material. As we have previously shown, the dynamics becomes particularly 
simple if we let the diffusive mobility ahead of the front be negligible % — » 0. In this 
case the dynamics ahead of the front, represented by Eq. (|2.4|l , is halted. This makes the 
A parameter unnecessary. This is not a particularly physical assumption, but it does only 
cause a small quantitative change in the pattern formation as compares to M = 1, as we 
have previously shown [4]. The existence of analytical solutions, however, makes M = 
an attractive choice. If we then assume that the front moves into material which has an 
equally mixed initial concentration, the morphology of domains formed depends only on 
the non-dimensional front speed (U). Interestingly, an analytical expression can be found 
for the non-dimensional domain wavelength (L) when U is small and constant: 

4(v / 6+61n(2- v / 273)-3 



The derivation of this law is given in [4]. Because it is convenient to work with this ana- 
lytical solution, we will focus on the condition where the initial mixed material contains 
equal parts of ^4-type and B-type material. 

To complete our model we require that the front moves with a time dependent speed 

U{T) = -j=, (2.8) 

appropriate for a diffusive velocity. The parameter C now becomes the only free param- 
eter in this model. This parameter will therefore determine the spacing and width of the 
resulting domains. 
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3 Derivation of Liesegang Laws 

We will now show that for certain parameter choices the model presented in the previous 
chapter will result in the formation of Liesegang patterns. We will do this by deriving 
the Liesegang laws of Equations (|1.2[) and jl.3} directly from the model. We will 

conclude this section by deriving a Matalon-Packter like analytical expression for the 
spacing coefficient. 

We first recognize that domain production at any point in space can only occur after 
the phase-separation front has passed. From our previous paper we observe that once an 
enslaved front has passed a point, domain growth or nucleation occurs very rapidly [4]. 
That is, domains form and grow near the front. The position of the front at time T is 
found by integrating the front-speed: 

X(T) = / U(T)dT= / — dT = 2CVf=aVf. (3.1) 
Jo Jo sjT 

From this we see that if domains are formed, then the nth domain will be formed at time 
T n at position X n which is proportional to \f%, in agreement with the time law of \\.\) . 
This is simply a result of the imposed front speed. 

Deriving the width and spacing laws is more interesting. As we found in our previous 
paper [4] and reproduced in Eq. (|2.7|), a front moving at constant speed produces domains 
of a predictable wavelength. In the process of deriving Eq. (J2.7|) , we discovered domain 
growth has two distinct stages: first is a formation stage where a domain nucleates at the 
front and grows as it is pulled along with the front until a new domain nucleates; second 
is an expansion stage where the just detached domain grows due to deposition of material 
excluded from the newly forming, opposite type domain. For a constant speed front 
moving into material which is initially symmetrically mixed and of negligible mobility, 
each stage accounts for growing half of the domain's final width. In essence, it takes two 
domain-type switching cycles to completely form a stable domain, and the domain width 
grows at half the front speed. 

The width of a domain as it detaches from a constant speed front is half of its final 
width, which is then a quarter of the constant speed domain wavelength for a front at 
that speed. If the front speed is changing, the width at detachment W dct is one quar- 
ter the domain wavelength predicted in Eq. (|2.7[) for a front moving with instantaneous 
speed U at detachment. This statement implies the assumption that our nucleation the- 
ory that was derived for constant U can also be applied for time-dependent U. This is a 
non-trivial assumption but it is justifiable be the excellent agreement of our theory with 
direct numerical simulations of the full PDE. Using Eq. (|2.8|) to replace the front speed 
dependence with time, then using Eq. (|3.1[) to replace the time dependence with position, 
we can predict the width W det of a domain as it detaches from a front when this front is 
at position X: 

W det = h(U) = = = ^X= 6X. (3.2) 

4 v ; 4U 4C 8C 2 H v ' 
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-- Final Domains' Width & Position 

Front- Attached Domain Width 
— Detachment Width W" 




Figure 1: This figure shows three important aspects of Liesegang pattern formation. The dashed line represents 
the distance from the previous domain in final morphology. We label the position of the first interface of a 
domain n as Xj|. Therefore the dashed line is X — Xj] where n is the domain that exists at position X. The 
dashed line has slope 1 because the W and X axes have the same scale. The dotted line represents the width 
W of the front-attached domain as a function of the front position X. The slope 1/2 is explained in the 
text and in our previous paper on en slav ed front phase-separation [4]. The solid line represents the maximum 
front-attached domain width of Eq. \3.2\ for the position X of the front. When a domain forming at the front 
reaches W det (X), it detaches from the front, and a new domain is nucleated at the front. The dotted, dashed, 
and solid lines all intersect at the points where domain type switching occurs at the front. Note that when a 
domain n becomes detached from the front its first interface becomes stationary at XjJ. 



Note that this already resembles the Liesegang width law, but does not relate the final 
domain width to the domain position. For that we must consider the growth of the 
domain after it detaches from the front. The information presented in this section thus 
far is graphically shown in Fig. [T] which represents the growth and final morphology of 
Liesegang pattern formation. 

To derive the width law, recall the formation and expansion stages of domain growth 
mentioned previously. A domain forming at the front increases its width at half the front 
speed. Due to mass conservation, the detached domain directly behind it increases its 
width at the same rate. The width of domains at detachment is predicted by Eq. \3.2\ . 
The final width of domain n is then its detachment width, plus the detachment width of 
the next domain of opposite material, which we count as n + 1/2. From this we get: 

W n = Wt et + W% 1/2 = flX d n et +pX d n % 2 . (3.3) 

The position of the first interface X° of domain n can be easily found from the detachment 
position X„ by subtracting the width W f f ef of the domain at detachment: 

X° = Xf t -Wf t = (l- j 6)Xf t . (3.4) 

As evident in Fig. [I] the position X° +1 , 2 of the first interface of the next domain, which 
will be of opposite composition, is simply the position of the current domain's second 
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interface, and is found by adding the current domain's width to its first interface position: 

X° +1/2 = X°+W„. (3.5) 
Combining Equations d3.3[) , (I3.4D , and H3.5D , we recover the Liesegang width law: 

We now derive the spacing law. As seen from Eq. (|3.5[) , the position of the first interface 
of the subsequent domain scales by a constant factor: 

X° +1/2 = X° + W„ = (1+ 7 )X°. (3.7) 

Two subsequent domains increase the index n by one, and their positions scale by that 
same factor squared: 

X° n+1 = (l+7)X° +1/2 = (l + 7« • (3.8) 
This gives the Liesegang spacing law 

X ° +1 = (l+V) 2 = - = I (3 9) 

The position of the band is taken as the boundary between the band and the previous 
gap. Due to experimental reasons, the location of a Liesegang band is often considered 
to be at the center of the band, i.e. X n = X° + W„/2. Our derivation of the Liesegang 
time, spacing, and width laws uses the fact that there is some freedom in measuring the 
positions of the domains X n , as long as all domains are measured consistently. It is clear 
that the scaling law holds whether one measures the position on the leading or trailing 
edge of the domain edge. Note that if one definition of the position given by 

X n = X° n +vW n , (3.10) 

obeys the Liesegang laws, then all these definitions of the position will still fulfill the 
Liesegang spacing laws because: 

Xn+l _ X° n+1 +vW n+1 _ (l+7v)X° n+1 _ X° +1 
X n X°+vW n (1+7V)X° X n ■ ^ ' 

From Eq. {3.9) . we determine that the Liesegang spacing coefficient for this model is: 

1 



(Y /4C 2 -iy 



1. (3.12) 



Note that this expression replaces the Matalon-Packter law, but bears no apparent resem- 
blance to this phenomenological law. For values of C > a/Yq / 2, the Liesegang spacing 
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coefficient p is greater than zero, and many domains will form, generating a Liesegang 
pattern; C cr = \/Yo / 2 sa 0. 1247 is a critical value. For C below this value no Liesegang pat- 
terns will be formed. Instead only two domains continue to grow and no nucleation of 
new domains occurs. A graphical representation of Eq. (|3.12|l is shown in Fig. [3] The au- 
thors are not aware of any previous purely analytical expressions that accurately predict 
the spacing coefficient for any Liesegang pattern producing models. 

We have now derived all of the laws for Liesegang patterns. Surprisingly there are 
no free parameters and our analytical expressions completely determine the patterns. It 
remains to show that the approximations made in this analytical derivation do not signif- 
icantly alter the results. To do this we compare our analytical results to direct numerical 
simulations of our model. 

4 Numerical Method and Results 

In our previous paper we present a one-dimensional lattice Boltzmann method (LBM) 
simulation of our model for enslaved phase-separation fronts moving with constant speed 
We did this by creating a lattice with spatially dependent control parameters, where the 
parameters change abruptly at the location of the front. We then advect the material 
at speed u across this stationary front, using carefully constructed inflow and outflow 
boundary conditions. A Galilean transformation x' — > x — ut of this simulation returns the 
original model. The simulation was designed this way, because a slow moving abrupt 
front on a discrete spatial lattice would be stationary for long periods followed by an 
instantaneous change of position. On the other hand, since the current of material in 
a LBM simulation is represented by continuous distribution functions, a constant drift 
speed u in the drift diffusion of material can be made arbitrarily small. For details on 
the development of the LBM simulations for enslaved fronts please refer to our previous 
paper [4]. Here we will only discuss the changes required to implement this method for 
fronts moving with a non-constant speed. 

4.1 Changes for Diffusive Fronts and Details of Implementation 

The method we presented was designed with constant speed fronts in mind, but there 
was no requirement made that the front speed, and therefore the material advection 
speed in the LBM simulation, be constant. To implement our model for the diffusive 
speed front of Eq. (|2.8|l we make two changes. First, the simulation is started at time 
T = To > to ensure that the advection speed is finite. Second, the advection speed is 
recalculated according to Eq. (|2.8|) at every iteration. To ensure the analytical solution 
of Eq. (|2.7[) is applicable, we must use front speeds that are sufficiently slow. The nu- 
merical verification of Eq. (|2.7[) is given in [4, Fig. 7] and shows that fronts moving at 
non-dimensional speed of U < 0.001 produce domain sizes that agree very well with the 
prediction. This being the case, we begin our simulation so that the initial front speed is 
Uq — 0.001, by setting our simulation start time Tq = (C/ Uq) 2 . 
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Figure 2: Final data from LBM simulation for C = 0.4. Subfigure (a) shows an example concentration profile 
of Liesegang pattern morphology at simulation completion. The offset is due to requiring an initial finite speed 
and dynamic growing of the simulation. Subfigure (b) shows the fitting the domains widths found in (a) to 
determine p(C). Some domains are ignored because they are not formed by the front, but are instead artifacts 
of the simulation procedure. The width W n of domain n is fitted to the equation W„ = Wg + Q(l + p)" to find 
an experimental value of p to compare to the theoretical prediction of Eq. J3.121 . For these data the fitting 
values are W = -1.92468, Q = 23.7999, and p = 0.212093. 



The simulation is initialized with random concentration fluctuations around the sym- 
metrically mixed concentration value of ( >, = 0. This results in spinodal decomposition 
in the region behind the front. These small domains serve to buffer the domains formed 
by the enslaved front from being adversely effected by any anomalies due to the outflow 
boundary condition. To increase simulation performance the number of lattice sites is 
initially rather small, and is dynamically grown to maintain a buffer of small domains. 

4.2 Measurements and Calculations 

Our goal with these simulations is to verify our analytical predictions of Equations (|3.6|) 
and (|3.9|l for the Liesegang laws. We accomplish this by performing simulations with 
different values of C, measuring the width W n of the nth domain formed, and performing 
a numerical fitting to the equation 

W n = W + Q{l + p)\ (4.1) 

where Wo, Q and p are the fitting parameters. This is an alternative form of Eq. (|1.2[) with 
a substitution provided by Eq. (|1.3|) . 

The domain widths are measured by interpolating the sub-lattice position of the zero- 
crossing of the concentration at an interface, and determining the distance to the sub- 
lattice position of the next interface. To ensure that there are sufficient domains to per- 
form meaningful fitting to Eq. (|4.1[) we track the number of switching events; after 19 
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Figure 3: Comparison of the theoretical prediction for p as a function of the free parameter C shown in Eq. t3.12t . 
This expression replaces the Matalon-Packter law for this model. The results show excellent agreement of the 
simulation results with the assumptions made in the derivation of the Liesegang scaling laws for this model. 
For C<C cr no new domains were nucleated. These data were taken from a series of simulations, the details of 
which are explained in the text. 



switching events have occurred the domain widths and other data, such as the concen- 
tration profile, are written out to disk. As explained in the previous section, the two 
domains directly behind the front are not completely formed, and are not used to find 
the experimental p value. The first domain formed by the front may also be ignored as 
it can sometimes be induced by the very strong dynamics of spinodal decomposition. 
This typically leaves 16 alternating domains of A and S-type material which can be used 
for fitting to Eq. (|4.1|) . Example simulation data and p fitting are shown in Fig. [2] for the 
C = 0.4 data point shown in Fig.[3l The concentration profile is shown in Fig. 2(a) and the 
width fitting is shown in Fig. |2(b)[ 



One additional note: the concentration profile of Fig. 2(a) seems to show a very sharp 
interface. The actual interface width covers approximately ten lattice sites, and is in 
agreement with the bulk stability requirements outlined by Wagner and Pooley in [11]. 
The use of the minimum stable interface width increases simulation spatial efficiency. 
Additionally, the use of smaller interface widths allows a larger time scaling factor s to be 
used which further increases simulation efficiency by calculating a larger time step each 
iteration. 

Numerous simulations of the type shown in Fig.|2]were performed for different values 
of C and the corresponding p values were measured and compared to the theoretical 
prediction of Eq. (|3.12|) . These results are shown in Fig.|3]and show excellent agreement. 
This remarkable agreement suggests that approximations made in the derivation of the 
analytical results have negligible effect on the final results. This is notable, as no other 
Liesegang pattern forming model has an accurate analytical prediction of the Liesegang 
scaling laws. 
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Figure 4: Example concentration profile showing Liesegang patterns formed by an enslaved front moving through 
material with an initial concentration <t> (fl = — 0.3, and a front speed scaling constant C = 0.5. The profile is 
shown in subfigure (a) at simulation completion. Subfigure (b) shows fitting the band („4-type) and gap (B-type) 
widths separately to find p(C) for initial material consisting of a majority B-type. Fitting parameters for the 
bands are W = -1.41736, Q=22.4956, and p_4=0.151002. Fitting parameters for the gaps are W = -0.44,9771, 
Q = 36.1552, and = 0.15435. This spacing coefficient is higher than for material which is initially equally 
mixed <t>j n = 0. As can be seen in Fig. [3] the spacing coefficient for C = 0.5 is p~0.12 in the case where <J>,„ = 0. 



5 Outlook 



In this paper we have shown that a model for an enslaved phase separation front which 
moves at diffusive speed U=C/ VT can be used to produce Liesegang patterns. We have 
done this by deriving the Liesegang time, width, and spacing laws from the model. Our 
analysis includes the determination of the Liesegang spacing coefficient p as a function of 
the front speed parameter C. In doing so we determined the values of C where Liesegang 
patterns are produced, and verified this with direct numerical simulations of the model. 

For this paper we have chosen the initial material composition to be symmetrically 
mixed. This corresponds to <£>i n = and allows us to use the analytical result in Eq. (|2.7|) 
to determine the Liesegang laws and their constants of proportionality. This model, how- 
ever, can generate Liesegang patterns for the range of initial concentrations between the 
spinodal concentrations. An example is shown in Fig.|4]of a Liesegang pattern for asym- 
metrically mixed initial material generated by this model. The parameters used in this 
simulation were C = 0.5 and 0,>, = —0.3. We expect that it will be possible to determine a 
spacing law for off-critical mixtures, but that is outside the scope of this paper. Addition- 
ally, the production of Liesegang patterns by enslaved fronts which move into material 
with non-zero mixed mobility should be considered. We expect the results to be qual- 
itatively similar, but we are not sure than this case will be amenable to an analytical 
treatment. 



13 



Acknowledgments 

We would like to think S. May and G. Kaehler for helpful discussions. E.F. acknowledges 
support by the National Science Foundation under grant DMR-0513393. 



14 



References 

[1] T. Antal, M. Droz, J. Magnin, and Z. Racz. Formation of liesegang patterns: A spinodal 
decomposition scenario. Phys. Rev. Lett., 83(15):2880-2883,Oct 1999. 

[2] T. Antal, M. Droz, J. Magnin, Z. Racz, and M. Zrinyi. Derivation of the matalon-packter law 
for liesegang patterns. /. Chem. Phys., 109:9479-9486, December 1998. 

[3] B. Chopard, P. Luthi, and M. Droz. Reaction-diffusion cellular-automata model for the for- 
mation of liesegang patterns. Phys. Rev. Lett., 72(9):1384-1387,FEB 28 1994. 

[4] E. M. Foard and A. J. Wagner. Enslaved phase-separation fronts in one-dimensional binary 
mixtures. Phys. Rev. E, 79(5):056710,2009. 

[5] F. Izsak and I. Lagzi. Simulation of liesegang pattern formation using a discrete stochastic 
model. Chem. Phys. Lett, 371(3-4):321-326, APR 7 2003. 

[6] K. Jablczynski. La formation rythmique des pecipites: Les anneaux de liesegang. Bull. Soc. 
Chim. Fr., 33:1592, 1923. 

[7] L. Jahnke and J. W. Kantelhardt. Comparison of models and lattice-gas simulations for 

liesegang patterns. Eur. Phys. }. Special Topics, 161:121-141, JUL 2008. 
[8] R. E. Liesegang. Ueber einige eigenschaften von gallerten. Naturwissenschaftliche Wochen- 

schrift, ll(30):353-362, 1896. 
[9] Z. Racz. Formation of liesegang patterns. Physica A, 274(l-2):50-59, DEC 1 1999. 
[10] K. H. Stern. The liesegang phenomenon. Chemical Reviews, 54:79-99, January 1954. 
[11] A. J. Wagner and C. M. Pooley Interface width and bulk stability: Requirements for the 
simulation of deeply quenched liquid-gas systems. Phys. Rev. E, 76(4):045702, 2007. 



